Stripes on a 6-Leg Hubbard Ladder 
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While DMRG calculations find stripes on doped n-leg t-J ladders, little is known about the possible 
formation of stripes on n-leg Hubbard ladders. Here we report results for a 7x6 Hubbard model 
with 4 holes. We find that a stripe forms for values of U/t ranging from 6 to 20. For U/t ~ 3 — 4, 
the system exhibits the domain wall feature of a stripe, but the hole density is very broadened. 

PACS numbers: PACS 
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The nature of the ground state of the two dimensional 
t-J model continues to be controversial, even if one re- 
stricts one's attention to numerical simulations. Density 
matrix renormalization group (DMRG) calculations l] 
for n-leg t-J ladders (with n ranging from 2 to 8) with 
J/t in the physical region of interest for the cuprates, find 
that stripes are formed when the ladders are doped with 
holes 0- In contrast, the Green's function Monte Carlo 
(GFMC) results of Sorella, et al., on square, periodic 
t-J lattices find a <i;;2_ a 2-wave superconducting ground 
state0. The aspect ratio of an n-leg lattice and the open 
boundary conditions used in the DMRG calculations may 
be sufficient to favor a striped state. Alternatively, the 
choice of the guiding trial wave function for the GFMC 
may bias the system towards a superconducting stateQ. 
It is clear that the t-J model is a delicately balanced 
system. Thus it is interesting to go back one step and 
study the underlying Hubbard model in which the onsite 
Coulomb interaction U is a parameter. 

Much less is known about the nature of the ground 
state of the 2D Hubbard model, which is computation- 
ally more demanding. Previous DMRG calculations 5] 
on 3-leg Hubbard ladders found that diagonal three-hole 
stripes with a linear filling density of unity formed when 
U/t was greater than 5. Early mean- field Hartree-Fock 
calculations for the 2D Hubbard model found that verti- 
cal stripes were favored for U/t less than of order 4 and 
diagonal stripes formed at larger values of U/t Q. These 
Hartree-Fock stripes have a linear filling density of one 
and their width is set by a coherence length £o ~ t/Ao 
with Ao the gap of the half-filled system. As U increases, 
£o decreases, approaching a lattice spacing at larger val- 
ues of U/t. Eventually, when U exceeds twice the band- 
width, the mean-field stripes disappear [3] . Note, how- 
ever, that the Hartree-Fock treatments do not include 
pairing and that their energy per hole is of order It too 
high. Here, we extend the DMRG study of stripe for- 
mation in the Hubbard model and in particular examine 
the question of whether stripes form on a 6-leg Hubbard 



ladder and how their structure depends upon U/t. Our 
results represent the first reliable ground state results for 
a Hubbard cluster large enough to exhibit an unambigu- 
ous stripe. From this single cluster we are not able to 
make conclusions about the infinite 2D lattice; on the 
other hand, our results strongly suggest that stripes are 
low-lying states in the Hubbard model, as they are in the 
t-J model. 

The cluster geometry was chosen with some care. 
Stripe-like behavior involving only two holes can be 
suggestive 8, 9], but may only represent a feature of hole 
pairs. One would like at least four holes in a stripe, which 
makes the system size beyond the current reach of ex- 
act diagonalization. Periodic boundary conditions, when 
the dimensions of the cluster are even, frustrate a sin- 
gle stripe0,0. Making one of the dimensions odd tends 
to force a single stripe. We choose a 7x6 cluster, with 
cylindrical boundary conditions, periodic in the y direc- 
tion and open in the x direction, and with four holes. Our 
previous calculations on the t-J model show that cylindri- 
cal L x 6 systems have ground states with four- hole stripes 
wrapped around the cylinder. The open boundaries nei- 
ther frustrate nor force the stripes, although if they form, 
the boundaries serve to pin them. One may worry that 
Friedel oscillations in the charge density, induced by the 
open boundaries, mimic stripes, but we find, for a number 
of reasons, that our stripes in the t-J model are inconsis- 
tent with the Friedel oscillation scenario |To). Sorella et 
al. compared GFMC and DMRG on a 6 x 6 cylindrical t-J 
cluster with 6 holes |3J; however, we expect 6 holes to split 
into a stripe and a pair, with resonance between them, 
obscuring any obvious signs of a stripe in the charge or 
spin densities. (A stripe having an odd number of holes 
would be frustrated in the y direction.) 

Here we begin by looking at the t-J model and then 
turn to the Hubbard model. A typical result for a 17x6 
t-J lattice with J/t = 0.35 and 12 holes, with cylindri- 
cal boundary conditions, is shown in Fig. 1. A maximum 
of rn = 5000 states were kept, for a discarded weight of 
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FIG. 1: (a) Charge and spin distribution on a doped 17 x 6, 
t-J lattice with 12 holes and J/t = 0.35. All the lattices 
that we will disucss have periodic boundary conditions in the 
vertical 6-site {/-direction and have open boundary conditions 
in the x-direction. Here, on the 17 x 6 lattice, three charge 
stripes have formed with 4 holes each separated by 7r-phase- 
shifted antiferromagnetic regions, (b) Same results, plotted 
differently. The open circles show the hole density (n^(Z)) and 
the squares show the spin density {S z (l))( — l) 1 . In addition, 
the stars show results for a 7x6 system with 4 holes, shifted 
to the center of the system. 

7 x 1CP 5 . One can see in Figure 1 that three charged 
stripes have formed and the spins form 7r-phase-shiftcd 
regions between the stripes. For an Lx6 ladder, the lin- 
ear charge density along a stripe is 2/3. Studies of longer 
stripes find that the preferred density on a long stripe 
is 1/2. Competing with the four-hole stripe, but higher 
in energy per hole, is a single pair, which also exhibits 
some stripe- like features. As seen in Fig. 1, there are 
two types of four-hole stripes: the stripes on the ends 
are (mostly) bond-centered, and the stripe in the mid- 
dle is site-centered. The energy per unit length of these 
two configurations is extremely close. The charge and 
spin density profiles along the x-direction are plotted in 
a more conventional fashion in Fig. 1(b). Since later we 
will be considering a 7x6 ladder for the Hubbard model, 
also shown is the local hole density on a 7x6, t-J ladder 
with 4 holes and J/t = 0.35. One can see that the open 
boundaries on the small system mimic the presence of 
the other stripes on the larger system. 

In Fig. 1, the presence of static magnetic moments 
makes the 7r-phase shift in the antiferromagnetic spin 
density as one crosses a stripe clearly visible. No ex- 
ternal magnetic fields were applied to the system, and 
these static magnetic moments are artifacts of the DMRG 
calculation, for which no spin symmetry has been used. 
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FIG. 2: Expectation value of the spin |(S z (i))l a t a site in 
the center of an undoped 17 x 6 t — J system, as a function 
of the error in the total energy per site. The error in the 
total energy per site is measured relative to a nearly fully 
converged DMRG calculation, accurate to about 10 -6 . Each 
symbol corresponds to a particular DMRG sweep, for which 
the number of states kept was m. Two sweeps were per- 
formed for each value of m. Although the exact ground state 
has S z (I) = 0, states with substantial static antiferromagnetic 
order can have energies only slightly above the ground state. 
The solid curve was generated by applying a staggered mag- 
netic field to the system, and represents the maximum value 
of | (Si* (0)1 possible for a given error in the energy (see text). 



However, we argue that the presence of these static mo- 
ments has very little effect on the validity of the re- 
sults. To see this, we show DMRG results for an undoped 
17x6 system in Fig. 2. The undoped system converges 
very rapidly with the number of states to, allowing us to 
study the convergence in detail. In particular, we con- 
sider (S z (l)) versus the error in energy for a particular 
number of states kept to. The exact ground state of 
this system, by the Lieb-Mattis theorem, is a spin sin- 
glet, which is spin-rotationally invariant, and therefore 
(S z (l)) = for all sites I. However, as measured by the 
spin-spin correlation function, this state has long-range 
antiferromagnetic order. The ground state is a superpo- 
sition of all the static antiferromagnetic states with all 
possible orientations of the order parameter. The reduc- 
tion in energy due to the macroscopic superposition of 
the states is very small. We see from Fig. 2 that there 
exist states with static moments of order 0.2 with en- 
ergies only 0.0002i per site above the ground state. In 
a doped system with some form of magnetic order, the 
behavior of the holes in each of the states with different 
spin orientations would be identical. It is reasonable to 
assume that the effect of the superposition on the hole 
behavior is slight. 

We can also analyze this effect in the undoped system 
by applying a staggered magnetic field. We assume the 
energy per site of the system, for small applied fields, 
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FIG. 3: (a) The ground state energy of a 7x6 Hubbard model 
with U/t = 12 and 4 holes versus the number of basis states 
kept in the DMRG calculation. The arrow indicates the ap- 
proximate point in which the stripe spontaneously forms, (b) 
The charge distribution (rih(l)) seen in the DMRG calcula- 
tion, labeled by the number of states kept per block m. 

varies as 

E(s,h) = E + as 2 - sh (1) 

where h is the magnitude of the applied staggered field in 
the z direction and s is the average magnitude of (S z (l)) 
in response. We minimize E(s, h) over s, keeping h fixed, 
to find s(h) = h/(2a), and hence E(s) = Eq — as 2 . From 
several very accurate DMRG simulations with various 
values of h, we find that for s < 0.1 — 0.2, this energy 
dependence is accurate, and for the 17x6 system a w 
0.0034. To describe the case where there is no applied 
field, and a finite value of s is considered an error, we 
define AE = as 2 , yielding s(AE) = (AE/a) 1/2 , which 
is shown as the solid line in Fig. 2. This should be 
considered to be an upper bound on the s(AE) obtained 
from DMRG, for which there are other sources of error 
in the wavefunction besides a finite value of s. 

We turn now to the Hubbard model. Our results 
for the Hubbard systems were performed using a new 
sine le-site" DMRG method[ll|, 

which performs better 
than the standard DMRG algorithm having two sites in 
the center when the number of states per site is more 
than two or three. Using this method, we have been able 




FIG. 4: The charge and spin distribution for the system of 
Fig. 3 



to keep up to 7500 states per block in some cases. Unfor- 
tunately, discarded weights are not informative and we 
do not report them. When performing DMRG calcula- 
tions on 2D clusters, one must deal with the possibility 
that the calculation will get stuck in a metastable state. 
For example, a striped state may be lower in energy than 
a state with two widely separated pairs, but the calcu- 
lation may only be able to tunnel between these config- 
urations when keeping very large numbers of states per 
block. In this case, one must repeat the calculations with 
constrained initial configurations, and compare final en- 
ergies. On the other hand, a calculation is particularly 
robust if one sees that it does tunnel between very dif- 
ferent states. Fig. 3(a) shows a plot of the ground state 
energy for a 7x6 lattice with U/t = 12 and 4 holes as a 
function of the number of basis states m. The hole den- 
sity distribution obtained at a number of sweeps, labeled 
by the value of m, is shown in Fig. 3(b). The initial con- 
figuration consisted of two separate pairs. As the number 
of basis states increases and the ground state energy con- 
verges, one clearly sees the stripe develop, with the "tun- 
neling" occuring for about m = 1200 states. From Fig. 4 
we see that, just as for the t-J model, there is a 7r-phase 
shift in the magnetization density across a stripe. One 
can see that as the stripe develops, the DMRG ground 
state energy decreases and that just as for the t-J model, 
doped holes on 6-leg Hubbard ladders can form striped 
ground state structures for U/t = 12. 

We have performed a limited study of the behavior 
of this system as a function of U/t. Fig. 5 shows the 
charge and spin densities for U/t ranging from 3 to 20. 
For U/t=8, we started the system with the holes as two 
separate pairs. We observed tunneling to the stripe state 
near m = 3600, considerably later than for U/t=12. We 
let this calculation continue until it used all the available 
memory (3 Gb), taking about 1 week of computer time, 
on an Athlon MP 1800+ processor. In this case, fur- 
ther sweeps reached a maximum of m = 7500, with little 
change in the charge and spin distributions. For U/t=6, 
we observed tunneling to the stripe state near m = 4000, 
again with a maximum of m = 7500. For U/t=4 and 
U/t=3, we did not try to observe the tunneling, instead 
starting with the four holes together in the center. Here, 
we found that a broadened striped configuration was sta- 
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FIG. 5: (a) The hole density (n h (l)) of a 7x6 Hubbard model 
as a function of the x coordinate £ x for various values of U. 
The solid lines refer to the list of U's on the right, in order 
by peak height. The dotted line refers to U = 20, which is 
out of order by peak height, (b) For the same set of systems 
the spin density (S z (l))(— 1) , showing the ir phase shift of the 
stripe. 



ble, up to the maximum number of states kept (m = 6000 
and 77i = 4000, respectively). For U/t=3, the charge dis- 
tribution is perhaps too broad to consider it a stripe, but 
the domain wall nature of the spin configurations config- 
urations is still fairly robust. However, since, as noted 
above, the spin expectation values are (useful) artifacts 
of the DMRG procedure, we do not attach much signif- 
icance to the specific magnitudes shown in Fig. 5(b). 
Fig. 5 also shows results for U/t=20. For this run, we 
started it in a state with the four holes together, but dur- 
ing the first several sweeps, keeping only m *=s 200, the 
holes partially split apart. Subsequently, near m — 1000, 
a definite stripe formed along with the antifcrromagnetic 
domain wall. A maximum of m = 6000 states were kept. 
For this large value of U/t, the stripe appears somewhat 
less stable than at U/t=12. We observe a very similar 
weakening in the stripe in the t-J model, as the value of 
J/t is reduced from 0.35 to 0.2 (the peak height in (n? l (/)) 



drops from 0.21 to 0.18). 

In conclusion, we have found that on a moderately 
sized Hubbard cluster, with cylindrical boundary condi- 
tions and doped with four holes, the ground state has a 
stripe. The stripe is narrow and well-defined for U/t=8- 
12. For smaller values of U, starting at U/t«6, the stripe 
broadens, until at U/t=3 the width of the stripe is the 
size of the system. At U/t=20, the stripe is somewhat 
broadened compared to U/t=12. The overall behavior is 
very similar to that seen in the t-J model. Although the 
open ends in one direction may encourage the stripe to 
form, the same roll may be played by neighboring stripes 
in larger systems. We find nearly the same hole density 
profile for a stripe in a 7x6 t-J system as in the central 
portion of a t-J 17x6 t-J system. For small U/t, the 
broad stripes we obtain are probably strongly influenced 
by the open boundaries. 
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Some of the reasons why our stripes are inconsistent with 
Friedel oscillations are: (1) the amplitude of the hole den- 
sity oscillations is very large, and does not decay signif- 
icantly away from the edges; (2) the hole density oscil- 
lations are largely independent of the length of an Lx6 
system; (3) four holes form a bound state, with a mea- 
sureable binding energy to breaking up into two pairs; 
(4) the amplitude of the hole density oscillations do not 
depend significantly on local potential terms on the open 
edges; (5) stripes form in the early sweeps of DMRG 
calculations in long systems, even when the hole den- 
sity near the edges is still precisely zero; and (6) suitable 
boundary conditions can induce longitudinal stripes. 
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